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The advantages of Goodman's Gamma, a measure of association, are discussed 
in i^eference to the Pearson coefficient of coi^tingency . Both theor ~ical 
and practical advanlfages and disadvantages are disqussed. » An empiri<^al com- 
parison the two measures shows tha^ gaitmia detects significant relation- 
ships which chi square does^ not, and that ga^T^.a is applicable to cases where 
cAi square is^ot. ^ \ - . . . ^ ; . 
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DON'T USE A CONTINGENCY COEFFICIENT, . USr) Gi^L^iJi ^ \. * 

' / ' ' ■ . 

In spite of the many warnings and difficulties (e.g., Lewis and Burke, 1949) 
involved in using chi square and. its related neasure of associat^n, the 
I earson 'continge^jcy coefficient, these statistics a^^ still used to-'t^fest 
hypotheses about order. In several widely used statistics textbooks, the 
examples used to illustrate a coritingertcy coefficient we;re two ordinal vari-' 
ables rather tJ^an two nominal variables. V7e think this 'practice is, mis- ^ . 
leading. Peri^ps usage of- the contingency coefficient stems from Walker and 
"Lev's {19S3^ p. 287) outdated statement that: ".Thi3 is no-t. a very satisfac- 
tory raeasure Of relationship, bUt under the circumstajices no better measure 
is available." However, tl\e situation has changed Since they wrote, 'and i't^ 
is the thesis of this paper that researchers should no Ipnger use the chi 
square and contingency*coef f icient for tests on ordinal variables. 

Goodman and Kruskal's (1954) gamma coefficient validly avoids the difficul- " 
ties that come with usir% the chi square and contin9ency^ coefficient on 
ordinal data. The gamma is a Measure of association for determining the^ 
relationship between twp ordinally scaled variables. Sophisticated computer 
programs, such as the Statistical Package for the Social Sciences ("SPSS) can 
compute this measure for contingency tables (Nie^ ^B^n^t, & Hull, 1970)r . 

Nevertheless, users of the SPSS or other packages, may be uneasy about using 
the gamma coefficient to; measure ordinal relationship^^ because they lack a 
convenient significance test. Although Goodman' and Kruskal (1963) have ' 
workfed out the sampling theory, the computation is complex. and has not been 
incorporated^ inrto standard computer program packages such as the SPSS. How- 
ever, we have developed a computer program to compute ^gamma, its standard 
error, its" 95% confidence inteirval, and therefore, its significance or non- . 
significance. We wiri be happy to send a documented 11 page user's manual, 
complete with FORTRAN 'listing,, to persons who write us. 

Now that the practical difficulty has beea solved, are there any theoretical 
reasons to prefer the gamma coefficient over the contingency coefficient? * 
The advantages s©em to lie with the gamma, because it has a simple and clear 
interpretation as a measure of monotonic association between two ordered 
variables. Gamma siinply reflects the percentage qf sample ^)airs which *are 
similarly ordered on two variables. More^precisely, if two persons are se- 
lected at random they are either tied on one of the two vari£ibles (x and y) , 
or one person has a higher scqre on both variables, or one person is higher 
on X, with tfie other being higher on y. ^Discounting tied ca^es, the gamma 
coefficient measures the difference between the percentage of similarly or- 
dered pairs. Or, it tells us the proportionate excess of concordant over 
discordant, pairs among all pairs which are fully discriminated or fully 
ranked. Furthermore, gamma has directionality varying from -1 to +1, and 
yqu don^t have to square a gamma t6 inte'rpret it.. 
^< 

In contrast, the contingency coefficient does not have a simple interpreta 
tion. Its square does not measure ^''explained variance, because in comparin 
tv/o ordinally scaled variables thjs concept> of variance is not meaningful. 
It has no sign; its upper limit varies with the number of rows and column 
in the contixigency tables; and it is not directly comparable to r, rho, 
or any othir cbrirelation coefficient. In fact, it is possible for one t 
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^ind a large chi square and contingency coefficient in situtations where the 
product-norent correlation is zero because of non-monotonic , non-lmear re- 
lationships. Tpis IS not what most researchers usually v/ant when they look 
for a^rnsjasure of association^. .Of course, if that is v^nat you want, by all 
means, use a contingency coefficient. Because of this insensitivity to or- 
dering, th'e contingency coefficient used as an inferential tool is less- 
powerful against population hypotheses of monotonia correlation than would 
be a test designed with such hypotheses specifically in" mind. Therefore, 
when some mono tonic correlation exists in .the population, the traditional 
con'tingency coefficient test is less lik'ely to reject the null hypothesis 
than, fc5r example, a test for the significance of gamma, which is ^designed 
specifically for such alternative hypotheses. 

v;$ -have recently compared the value of gamma and the coHteingency coeff i- 
.cient as^ methods' of describing the association between measures of leader- 
slnp at the U.S. Miritary Academy and four criteria of post-Academy officer 
performance (Butler, 1973) . The study wafe based on the Classes of 1961 
through 1965. In all, 112 contingency tables,- either 3x5 or 2x5, were de- - 
veloped. For 34 of the'se tables, it was impossible to compute a chi square 
cind 'the contingency coefficients because the expected frequencies were too 
smair*. In 18 others, it was necessary to combine columns and rows. In con- 
trast to'this, the* cell N's were togt sm^ll to evaluate the significance of 
gamma in onl^ 6 of the tables. This is so because of the sampling theory 
developed by Goodman and' Kruskal <1963) , which allows a researcher to eval- 
uate the significance of a gamma for a wide range of sample sizes. Goodman 
and Kruskal developed 'large sample theory for the gamma coefficient, and its 
standard errors and significance tests are based upon the findings that 
gamma has an approximately nomal sampling distribution for a sufficiently 
large sample size.^ From sampling experiments reported by them and 'by 
Rosenthal (1966) , it is apparent that large sample statistics can be safely 
applied when there are an average of four cases fot every cell in the teible. 
For example, in a 3x4 contingency table, a sample size of 50 is adequate. 
Goodman and Kruskal (1963) also give a "cons.ervative" standard error which 
can be used when the sample ''size is smaller than an average of four cases 
per cell. When the average frequency per cell is between two and four, the 
"conservative" standard error apparently produces acceptable results. Our 
computer program computes tfVie appropriate standard error, (depending on *sam- 
ple size. In our study, when both statistics were Computable, nearly 70% of 
the comparisons showed that gamma was larger thah the contingency coeffi- 
cient, often by a .fairly sizable amount. There were seven cases where the 
gamma coefficient was significant at the 5% ievel, whereas -the contingency^ 
coefficient was not significant. As discussed earlier, the probable reason 
fo]Q this is that gamma is more sensitive to order. Furthermore, the gamma 
allowed a clearer interpretation, indicating the size of the monotonic asso- 
ciation between two ordered varicibles. However, a significant contingency 
coefficient showed that there was a significant association among^the vari- 
cibles, but not whether the relationship was monotonic, curvilinear, hyper- 
bolical, etc. Further tes£s would have to be made to determine the type of 
relationship v;hen using the contingency coefficient. 

Tv/o 9ther recen,t studies* have compared the sampling stability of the contin- 
gency coefficient with other measures. Whitney (1972) compared the chi 
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square statistic and a rank order coefficient developed by Kendall. He 
showed that the rank order test is more pov-erful, proviqled that the under- 
lying relat4.onship between the two variables is linear or monotonic. 
Sarndal (1974) evaluated fifteen measures of association and concluded that 
the coeff icie?it of contingency was among the most biased (pp. 185-186) . In 
conjunction v/ith data from the present research, one can conclude that the 
contingency coefficient is never the coefficient of choice. 
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MEMORANDUM FOR: 



RESEARCHERS AND PROGRAMMERS ' . 

OFFICE Of the DIRECT05L OF INSTITUTIONAL RESEABlCH 



SUBJECT: ' Program GAMMA 



1. Purpose : To compute GAMMA, a measure of association. 

2. Background ; ' , ^ 

a. Goodman and Kruskal (1963) proposed GAi^MA, measure of association', 
whicl^ reflects the degree of order in a contingency table. Unlike the' 
product mojnent correlation coefficient, which is based on the principle 
of least squares the G7VMMA coefficient optimizes the prediction pf ordo , 
The logic of tide GAJ^IMA coefficient is as follows: if two persons are 
selected at random, they are either ti^ed on one of the two variables (x and 
y) , or one person has a ^higher score on both variables. Or one person is 
higher on x with the other higher on y. Discounting tied cases', the two 
variables are positively associated if there is a high probability for 
variables x "^nd y to rank order two persons in similar ways. The two 
variables have a negative relation if variable x; tends to rank-order 
persons opposite to the way variable y rank-orders them. The GAMMA, 
coefficient is b^sed on a direct count of the aumber of similcurly rank- 
ordered pairs (S) minus the number bf dissimilarly rank-ordered pairs (D) . , 



GAMMA = 



b. The GAMI-IA coefficient has a rationale wliich is similar to the 
rationale for the »tau coefficient, prbposed by Kendc^ll (1970). The t*au " 
coefficient is influenced hy tied scores^ whereas GAMMA is not. GAMMA 
is a parjd.cularly valuable measure to use, therefore, when there are 
many tied^cores, such as one finds , in contingency tables. 

c. Goodman and Kruskal (1963) developed large sample theory for the 
GAMMA coefficient. Standard errors and significance tests are based upon 
the fact that GAMMA has an approximately normal sampling 4istribution for 
a sufficiently large sample size. From sampling experiments reported by 
them and by Rosenthal (1966) , it is apparent that large sample statistics 
can be safely applied when there are cin average of 4 cases for ev^ry 
cell in the table. For example, in a 3x4 co^ntingency t'aiile, a sample 
size of 50 is' adequate. For a 5x5 table, N should be ^t least 100. 
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SUBJECT: Prog^ram GAM2-L?V 

d* Goodman and Kruskal (1?63) also give a "conservative" standard 
eri^or which can be used when the sample size is small. When the average 
frequency per cell is between 2 and 4, the "conservative" standard error 
apparently produces acceptable results. ' 

3. Description of the Program ; * ' 

Inputs: The pr'ogram reads contingency tables (up to* 20x20) • 

Outputs: and .(1) prints each with a label 

(2) computes GAi4iMA 

(3) if the sample size is large enough computes' 
the standard' e^ror and 

(4) a ^95% confidence- interval for the population GAMMA. 

^ The 95% confidence interval simultaneously tests a number of hypotheses 
^ about -GAMMA. If the confidence interval does not include 0.00, then you 
can reject, the usual null hypothesis that GA>U^IA - 0,00, In other words 
when the confidence interval includes only positive non-zero values, you' 
have a "signif itant" positive as^ciation. At the same time, if the 
confidence interval goes from .7 to .2, you may reject the hypotheses 
that the two variables co:EC2:elate more than .7. 

4* How to Use the Program : First, sort all the tables on the same size 
(same number of roWs and columns) into one contiguous group-. 

Table Size Card: * ^ 

Col 1-3 NR - Punch t:he number of rows in the contingency table 
right justify.* 

Col 4-6 NC - Punch the number of calumns in the contingency table^ 
right justify.* 

Col 7-9 NT - If there are a series 6f tables 'with the same number 
^ of rows and colximns, punch the' number of such tables. 

If you oniy have one RxC table, leave Col 7-9 blank. 

Contingency Tables: Each 'table will require 2 or more qards, .Each 
card lias the same basic format. 



*Right Justify: TT you hav^a 3 digit number, leave two b^ank spaces, 
then punch the 3 digit intoVhe right-most part qf the 5 4^git field. 
In general, if you have a shdrt number to go into a large space, move 
the number over to the right hide of the field, leaving 'the left side 
of the field blank. \ 



The table should be arranged «so that (reeling left to right) each 
successive column represents successively higner values of variable x, and 
each successive row (moving down from the top row of the table) represents '"^J'^H^ 
successively higher values of variable y. ThuS/ the first card you 
punch begins with the lowest ranked coll in the' entire table; the last 
card you punch ends with the highest ranked cell^^ In the data deck, place 
the cards in the order punched, with top row (lowest y-vaiues) first. For 
example, in the^f ollowing tcible: * * ■ • 





- Low 




Med 


Ki 








Low 


10 


N 

5' 


0 


1st 
' 2nd 


Card 




Med ' 


• 5 




11 ■ 


4 


Card 


• 


Hi 


0 




3 


12 


3rd 


Card 





There are 10 scores with low scores on both X and Y; 12 scores with high 
scores on both a and Y; and so on. Start each new card with a label, and 
punch the call frequencies within one row of the contingency table, 
allowing five digits for each call frequency. Start each new row with 
a new card. . . * • , 



Col 1-5 

COL 6-10 
COL 11-15 
COL 16-80 



FORMAT FOR EACH ROW OF EACH TABLE 



Punch a 5 character label into each card. This label 
ider^tifies the table for you in the ourput, 

.Punch Row i. Col 1 frequency, right justify 

« 

Punch Row i. Col 2 frequency, right justify 
Punch Row 1, Col 15, frequency, right justify 
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Symbolic Deck order: 



Control Card 1 



Table 1 



Example: 

,2 2x2 table, with 
2 cards per table. 



Tablfe 2 



Last 2x2 table 



Control Card 2 



Table 1 of the 3x5 tables. 



2nd 3x5 table 



. etc*, more 3x5. table 



0001 $ SMJMB 

0002 $ IDENT 



EXAMPLE OF JOB SET-UP 
(A) Control Cards Which Preceed Program 

XZ971,40 

,MAOR/DL , ' 



0003 $ LIMITS 02,25K,,1000 

0004 $ OPTION • FORTRAN, MM> 



0005 $ FO^TY 
0006$ INCODE 



USMA 



i 



GAMMA Program 
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,-0"8-73'' 15.417 '(^OOnMAN-KRUSKALS GAMMA _ > _ 

C GOODMaN-KRU'sKAUS GAM.^^^ • . . 

r ' rODEO BY PRiEST-MAnR OCTOHER 1973 , „ ^ , . 

r " THIS PROGRAM PRODUCES A MEASURE OF* ORDER ASSOC I AT I ONCALLED GaMkA 
r ThE GAMMA COEFFICIENT WAS PROPOSED PY GOODMAN. AND KRUSKAt „ , - 

MEASURES THE CONDITIONAL PROBABILITY THAT TWO RANDOMLY CHoSENj^'SRSO NS 
^ JIlL RAS'K IN THE SAME wAy ON VARIABLES X AND y. MI^US THE PROBABILITY. .. 
•C THAT THEY WILL RANK IN OPP-OSITE wAyS*'ON VARIABL X AND VARIaSL Y j 
C GIVEN THAT THEy ARE NOt TIED ON EITHER VARIABLE, 

C THE PROGRAM COMPUTES THE GAMMA COEFFICIENT . ... 

C_ ITS STANDARD ERROR ' , r 

C AND. A :9t> RERCENT CONflDENCE IMTERVALAL FOR ThE | 

r POPULATION GAmMmA „ • — 

C • IF THE SAMPLE I? LARGE ENOUGH THE BEST ASYMPToTiC 

r StAMDARDAR ERROR WILL BE 0 COMPUTED '- - -■ 

Cc' IF THE SAMPLE SlZE IS SHALL? THEN : 

/4 - THP pROGRAM COMPUTES AuCO^^ISErvAT I VE» StAnDArD ERRORi — . 

c '■ It IS THE UPPER Bound For the more precise standard 'ERRy 

r ' FORmULA ' ' . . ^ 

C~~LF THE SAMPLE SI^E I.S foO SMALL, THE PROGR.AM _PR I NTS_.A .MESsAGETO-THAT-^eFFE 

C -CT ■ . -: I 

C HOW TO USE THIS PRoGRAH ' :" ' 

C FIRST OF aL'l; there IS ONE CONTROL CARD FOR EACH TABLE ' ' 1 

C~OR SERIES OF tIbLES, THE CONTROL CARD TELLS T^HE PROGRAM JHE nUmBER ; OF j 

C ROWS, THE NUMBER oF COLLUMNS IN THE NExT TABLE OR SERIES Of TABLES ■ 

C THE FIRST CONTROL CARD ALSO TELLS tHE ProGRAMTHE NUMBER OFJABUeS HA'Vc JH 
a _ SAME DIMENSION. tHES6S THREE PIECES OF I NFORMAT I ON ( P A.R AmcTErS > ARE-. CALL 
C NR, nC, NT By the pROGRAM. PUNCH THEM INTO A SINGLE'-CARD : 

C ' AND ALLOW 3 DIGITS FoR EACH NUMBER P IGHT. •JUST I Fy. EACH .NUMbEr .PUN^CHED', 

C I 

C ALWAYS PUNCH A FivE DIGIT IDENTIFICATION NUMBER AT THE BeGINNInG Ofi EACH 
C— CARDi THIS FIVE DlGlT I D.ENT I F I C At I ON NUMBER WILL BE USED TO HcLP YOlU „J 
C IDENTIFY THE PARTICULAR TABLE WHiCH PRODUCED THE GAMMA ^- ■ • , ^ ! ! 
C THEN CONTINUE PUNCHING EACH CARD WIIHTHE CELL FREOUENCIES ONelROW "0^ _ 
C THE TABLE. ALLOW ^IVE DIGITS. FOR E^CH CELL FREQUENCY, ANd RIGHt JUS.;TIFY' 

c_£ACH Number if there are fewer digits in- the cell frequency., ..^ — 

C " START EACH NEW ROW OF THE TABLE WITH A NUEWCARD OR SET OF CArDS' ,\ 
C EACH TABLE SHOULP BE paNCHED..AS ..A . SER I Es .QF .5,. D I G ll. NOMBErS 

C ■ ■ ^- . 

C : 



DIMENSION X(20*20) 
DIMENSION S(20f20)iD{20,20j, IrEG(4,4),y<^J 

12 CONTINUE 

READ(li,iOO,END = 99)__fclRV.Nc,NT_ 

1F<NT.EO;0)NT = 1' - ^ 
DO_4.0l._NtABLS5l|NT ^ 



C ' ' . V 
Z\ NR B. NO OF row's; .Nc s. NO. COLUMNS V 

C 

Do 400 I cIx'nR ♦ - ,. 

READhl»101,ENP = 99)A0><X(IiJ),JBl,NC ) 

O i_PDPMAT{A3.iOF5!0,i -• 

ERLC-T-T^-^ r-^f^r--^ 

^' ' ' ' ^ PROGRAM LIST {!) 



J 



J 
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'idODMA^.-ti-fUJSKAUS GA,<^1A. 

^lw?Li RA^Ik I,' THE SAME WAY ON VARIABLES X AND y.'MINUS THE pRQ^BlLlTY 
•^•THAT THEY -^lUL- RAV.'K IN OppOSIT^ WAyS ON VA-RI^rfL / AND VARIaBL ' , ^ 

IgIVEN THAT THEY ARE N'O-T T I E.t)-- 0-N EITHER VARI-BL'E.. '■ 
\ ^-Hc 'PROGRAM COMPUTES ^THE .GAMmA COEFFICIENT , ' . : ^ ■ . . 

>-f-^'^ us STANDARD ERROR • . . v c • 

♦ ' • . A'nD a ^'p" PERCENT CONFIDENCE INTERVALAU FOR IHE -j^- 

'pOPuUTlON' GAMMMA „ ^ --■^ — • 

iF^tHE SAMPLE jSfLATCE ENOUGH THE BEST ASYHPToTI-C ■ 

St andardar Error will bE*o ceMPpTSo v. - - ■ ■■ ■ - ■', 

C'^f'"' • ' ll THE SAMPLE SIZE ' I & SHALL'- T.HEN 



THE pROGRAH COHPUTSS AnCONScRVAT I vE" StAnDARD ERROr'J 

THh P^oGKg"^;;^ ^j^^pgp Bound ,foR the-hoRE precise standard frr^? 



^JFTHE SAMPLE SlZE'is\oO SMALL.^ THE PROGRAM PR I NTS ^A _MES:SAGET0_7HAT— EFFE 

' Vtp^T nr 'Ai'i 1 twErp IS ONE CONTROL CARD ^FOR EACH TABLE - ^ 

"orsEM s' *tI3CeL tSe CONTrSl'cARD tells THE PROGR-AM THE nUhBER OF 
' ROWS, tH^ NOHBER-yi^.COLLUMNS IN THE nE>cT TABLE OR SER I ES OF - TAbLES ■ 
'^t-bE F RSt"cONTROU^ARD ALSO TELLS, THE ProGRAMTHE NUMBER 0?' IT^Q^.^S H A ' Vc TH 
Ia-ME DIMENS'ION. THESE'S THREE PlECES^OF 1 NF ORH^T I ON ( PARAmEtErS > ARE. fcALL 
"T" WD KiC, mt Rv TwE PROGRAM. ■ PUNCJH THEM INTO A SINGLE"-CaRD v _ v . 

^aSdAlKJdIGUs Cor EAcS NUHBER.: R I GHT- -JUST I Fy EACH.NuMbEr PUNjrCHED". 



■'"i^'acSa?s punch a five Di.GIT I DENTI F Tc AT I On' NUMBER f THE BeGINnI nG "ofT EACH 
'•!^CARD. TtlfS FIVE D i 6 1 T " 1 DENJ I F I C At I ON Nui^ER wTllIBE USED TO, HElP YOi U 
" IDENTIFY- THc PARncuUR' TABL€ WHiCH PRODUCED THE OAMMA ■ ^ ! _ ■ 
r {hen CO^tISqE PUNCHING, EACH CARD WITHTHe CELL FREQUENCIES INONe.ROW OF 
r-m TA6L£- ALLOW ?^r^E^D I G I TS FOR tA-CH CELL . FREQUENCY , AND' R I GHt JUS: ? IFY 
c'^ACH NTiMBER iF THERE' tpt FEWER DIGITS IN THE CELL FREQUENCyn .4- 
r7^4xAR?.tA^C.^ N€^ ROwsOF THE TABLE W I TH ;A NUEWCARP OR SET'OF CArDj?. • 
l__:iACH" TA&lt S>^OUL.■D„BE PUNCHED. AS A. SERIES OF. 5 . D I G ll. NUMBErS ^ | ^ 

! -.•. '' ■•■^^N-^ ' \ ; " ^ — — I- 

■ - . DLMEwSii >:'(20--*2o')- ' , ^ " ' ' ' i 

• •^_^I>IM^Ns;o'■N^.»S{20/^> VDt^0'v2:<iinREGi-'^''^i'V(^J — . \r — 

1^-* CbNT-I'N'UE ' ' ■• k • „ ■ • ' . I . 

•!:2ll^READ.{.ll^iOO,END^9). tJRyNCjNT 



.:4r<NTV^o".o)NT=; 

_ nOfL AO l.:,NtABLS4«N T -■ e ^T^^ 



I 

1~ 



•._?L_f*f^tf. N_p.OF 'ro(^'s, ':H 9>o c.olumVs 

Do 4(ll0- I«1VNR ' . , • • - - -V ' 

• ' ^ READ{"li,i01,ENP = 9^)^C><XM»J>iJ»l|NC \ 

^ ^i JM M 3 » :toF\5iOv ^ 

•V f 1 V " •'* ■ ■ . • - ■ ' PROGRAM LIST' (1) 



A 



♦ ERIC • . ' [-*^ ^ . 
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11-08^73* 15,417 ^OOt^MAN-KRuSKALS* GAMMft 
, ' ^ ■ , ■ ' ■ . • 

PRINT i07,AD.VUri » J) 's) = l»NC) 
'l07 FO^«AT(iH' iA6,iOr6|l) 
•400 CQNt'InUE ' ' y 

DO'^S^ 1=1, • 
:.. .DO 33 J=1,NC ^ 
33 TX=TX*X( I , J) 
RC = N'9«NC 



'RC = N'9«NC . 

ireg is k nEvicE For controlling operations about a. pIvot in 

'- lREG{l.iU=.— 1-— ^ ~- ^ ^ 



A, MATRl; 



nREG<l»3) 
- lREG{2ai 
lREGt3>3jf 
. IrEG(2*4-7 
lR£G(3/*2y 
-IrEG(4'-27 
lREGt'4.4if 



— .^104)- F0(ViAT(3.j3) 
C • » 

_c: srARt ■- 

C ^ - ■ 





= nC 

=_nr — 

NC 



- DO .1 I=i,NR- - _ 
IrEG(1»2> = , I?l 

-lREG(2»'2y = I"!-- 
lREG(3».lif = 1*1 
-IREGMiIT -= -1*1- 

DO 2 jsljNC 
-iREGdM) ■= -J"J-1- 
, IREG(2'3) = ^♦l 
-JREg(3.4) = J-1.. 
IrEG(4,3) -Jtl- 

— Do. 4 K3l*4 



I A. = IREG{K,1) — ^ — 

IB = IREG(K#2) 
JA s IR£G(K*3) 

Jb = IREG('<,4) 
—iFdA.GT'. IB) GO- -f 0-4. 

IF(ja,gt;j8) gO to 

DO 5 i 1 = 1 A' IB 

"qO 5 JJaJAiJB 
.5'. Y(K)aY(KV + Xni'/jJ)- 
4 CoMt'INUE 

S( I, J).= Y(l)+Y(<*->' ^ 

D( I' J)=y/2)*Y(?) ■ 

.2-CO'MtInUE - 

1- CONTINUE - 

PS=0 

PD = 0 

P-SSi;0, — ^ ^ 



•A 



1 



1 



(2).v 



ERIC 



11 



14 



11-08-73' l5,4i7 ^OOOmAN-KRuSKA'LS GAMmA" ^ ^ 

PD3 = 0, ^ " ' •" • ' " ■ * 

• ■ ?SD = 0, ^ . . 

Do 32 Isl,N'R 
DOv 32 J=l i NC 

C PS IS tWg NUMBER Op PAl^^S Of PERSON S WHO A^E FfANK£D""S I H I LaJ^L Y~0.m 

- C-PD !S THr NUMBER OF D I SS I M I L AR . f^A I RS ON BOTH X AND Y VARlAjES 

ps=p«i*X( I » J)<»s( i Vj), ' ' ■ • ■• ' 
• . ?d=pd^x( 1 , J.)*D(r;j) . . ■ . 

■ PSS = PSS> S<l/j>*i^{I» J)»SM,J) • 

?DO=PDD* Ddt^^yil ,J)<^Dil,J)^ . . 

. ?SD = PsD*S(I»J)*^<"liJ)*D(rrj> 

32 CONTINUE ■ ■ > . . ' ^ ^ 

GA'^MA?(ps-PDi /^PD *PS) 

C.w^.l.AR3c ■ SAMPLH ASyHpjOTlC STANDARD ERROR ,_ 

ZEsDapS ♦T'DDepS -2 , 0*pS'»PD*'PSD ' ♦PD»PS9»PD 

: ZEETeS^RTrZEED) • . . . , 

'2E£ = (PD+/PS)*»2/'t4,D»ZEET5 . . , > 

'^..-Ifi (T^/RC-i .Ci:4j0) GO 70 34 „ -.Ijl'^^ 

' PRINT 35 ■ • . 

^35 roRhAT( '■' SAMPLE SJZE PER CELL IS NO.T LARGE E.V'frUGH FOR »» 

X " HOST POWERFUL ASYMPTOTIC TEST." ) 
Ifl (Ty/RC)-;Lr»2,0) GO -'TO-. 37 



C A MORE COMSERVATIVE STaJ^DARD^ ERROR FOR SMALLER S-AMPLES 

Z£EDk(PS*PD)/(2, OftTx»(liO-GAMMA##2 .j )'....,_.'. 

ZEEs ' SORtfTEED) 

34 CON'TlNUe „ 

• ErR=1.0/7eE 

-PRINT 105, GAMMA, ERR - .• 




105 FoRmATC 'i sample GAMMA, b «,F8,4,if STANDARD ERROR » i F8 , j 

— OL - sgamma*1i96<»err . „: : 

Wl- sGiiMMA-l,96*H9R 
PRINT 106,Ul7wL 



.106.FoRHAT( 'I PROBABL'YC ,V5) POPUL-ATION GAmMA" IS' BETWEEN n 

_X?'8t4i. V AND U.F8.4 ) ^ ■_ 

GO TO 401' 

37lPRlNT 3^.- . J). L . _.. , 

38 F0RMAT(> THE SAMPl^E SIZE PER CeLL' IS; TOO SMALL'O 
._PRlNT 402»GAmhA .* 

mo2-forhat(ih "/"gaMMa Is ir, F10,5> 

_4 01_CONTfNUE - [ , 

Go TO. 12 ' 

99 STOP ." ^ __ . ■ . ' ' 

' - End- . . •. 
i . 

!E7no' ' DIAGNOSTICS "I-N "ABOV'S''cOMpyL"AT I on T*!? 
)S_HERE:.USED-F-0R_THlS_C0MP.lL6XIpN__^k " ■ 




PHYSICAL PICTURE OF COXTRO;: 
• AND DATA CARDS FOR 
SAMPLE PROBLEMS 



CARDS 



TABLE SIZE 
CARDS ' 



/ 



//L_r_ :plC3. . ,0. U [f-l 



JL ^\l j-l-L-LJ-i-JU.'- - j^^-jjim l. I ,U.nni!M Vr,TT.ririfi.'.rn7ikir^W3 , 



J ; 3 i 5 fr 7 ,,-: ^ r. 1 « 0 ' J 

TTo^Q Q pia Q 0*0^3 P^^^t t Pi i L.I L. a— IJ -Jt 

1 1 ill ilu_tU w'kirj^M 



•Third 
Problem 



Second 
Problem 



/ 



^\ ^1 I _ ! I II ' f !.! I -J) - ^1-1 i 

! pW' ti?L I P t 1 1 ^ L-J !-J 1' 

; IXC K > ■ lU nirnkiruhTO ^ . 

y ( Jr[^H > i t !»a!ii w iW » »Iw u ?7 73l;« s »t?r^^ ^ 
2^ * • I 1 1'* 



First 
Problem 



USMA 
tONTROL 
CARDS 



» ' > ^ , .^IHpODE USMA , , t il I f J 

Oj^f g j ' , a } 4 r T7^*^Tir?rn~»]"»" » » » »» » » " ^* » " " ^' " " "-j j-'j-".?! 
Ir* / . ^^"i'h . V'*^''^^^Pt !i I L_LL_I_I— mJ 



^iHLxw^ . . r - , a > 4 ^ t > » » » n n n M t> IS 7t ;i » n » ;s ?s 71 3i » » » » » » » I LiUi 

^ ' . . . rrr- -. . t — JT ij «t 'It II M 1« la II 42^ 



4* 



f ' 

-J 



,0 I 

oTo'3 0 0 0 0 0 oM 0 o'o 0^ 



333 

M4 



222 



1 \\ 



2,22222 



86B 



11 VI 



3^3 

4444444 

5 5 5 5 5 5 S5^5 



3133 



S 8 8-3,3 8 
9 9 sjo 9 9 



6 8 &i 5 6* JO 
ONE 



1 I 1 
222 
3 3 
444 
55S 
8G* 



1 \ \ 
222 
33 
4 4>4 
355 



:3 0 0 

1 1 1 

2 2^2 

3:k 

444 



?T7Trrr7[n7imp min ? 



8 8 8 8 



999 



o'loo 
1 1 1 

222 
33 3 
444 



6 65668? 
TWO 



9S9 



880 



888 



9999999 



It U i) it il II 17 li 



0 00 

\ \ \ 



2 2 2 2 2 2 2 2 2 2 2 2 2 



555555 

UP 



0 0 0 0 

1 \ 1 



333 



555 

8G6 



333 



4444444444 



S55 
888 



THREE 



\t ?7 2q__/» i*< /* ^' — " ' 

0T?0l O'O 0 plO'O 0 0!0 0 0 0 C'Q 0 0 0 0 



1 1 1 



3333333^33 



71 u n 



888 



9 9 9 9 3,9 9 9 9 D 9 9 9 

^ >^ M «A M «■ «^ 



21 « jc1j7 n 71 



888 



in 



1 1 1 
222 



4444444 



5 5 5 5 5.5 5 5 5 



1 1 ijni 
2221222 



1 1 1 



1 1 1 



2222222 



3 3 3'3 3 3 3 3 3 3 3 3 3 3 3 3 



8 8 8 8 8 Gl8 
FOUR 



i5S5 
; 6JJ6 



444 



444 4 444 



444444 
555555 



r 



88 



5 5 5 5 5 5 5 5 5 5 

6 8 8!6 6 6 6 8 6iS 6 8 8|6 8_8, 

3 8 6'8 8 8 daii 



000 
1 1 1 
222 



9 9999 9 0 9 9 919 9 9 9 9 009 SS;299 

(37 « »,<J»«I U O Jl4 i> « «|»l5i M ii;>4 « W 



ERLC 



lb 
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c. Exairiple of Program P'utput 



> 



o o o o o o- 



o 



o o 



ErIci 



> t 











o 




o 




o 


n- 


Di 






tr 








Q 
























rt 






















»— ^ 








e> 










0) 




0) 




o 



























0) 

rt 

0) 

n- 

& 

CO 



I 1 



1/ 



^ >• m 
o 3: o 

r" 

•<<o 

> o 
n 

0 i 

no 

JX> -s) ^ 

o 
o 

^ >- 
3: :2 -c^ 
^3: a - 
o 

■ >o 

CO 

1 m 

COX) H»- 
— 1 o - 
m 

m ! 



ooo3:wro*-*o3:x3>> 



^ > 

o : 



I 



o coUJ>oocDvn>'OVJio 
It • o n : o o 



o 



o n 



o 

I no 
! cr o 

I r-* r- - 

. ^-^h*'•---^o^--- 
' o ro. o ^ *^ 

-u w :s h*- o*;^ 12: o ui o 

i i ' I XI . 

I CO ; CO 

1 m ! I m I 

I CO XI 1 'COX) I 



I 



I 



o ' 



o a> ! 



o 
ro 

i 



m 

m ■ 



o. 

C33 

I 



m j 

z 

I 



I 



5>l 

t3 



I 



I O 
O C^ 

r 



o 

V3 



00 
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